Non-singular cloaks allow mimesis 



Andre Diatta u and Sebastien Guenneau 2 
September 6, 2010 

x * University of Liverpool. Department of Mathematical Sciences, 
M.O. Building, Peach Street, Liverpool L69 3BX, UK 
Email address: adiatta@liv.ac.uk; 
2 Institut Fresnel-CNRS (UMR 6133), University of Aix-Marseille, 

case 162, F13397 Marseille Cedex 20, France. 
Email address: sebastien.guenneau@fresnel.fr; guenneau@liv.ac.uk 

Abstract 

We design non-singular cloaks enabling objects to scatter waves like objects with 
smaller size and very different shapes. We consider the Schrodinger equation which 
is valid e.g. in the contexts of geometrical and quantum optics. More precisely, 
we introduce a generalized non-singular transformation for star domains, and nu- 
merically demonstrate that an object of nearly any given shape surrounded by a 
given cloak scatters waves in exactly the same way as a smaller object of another 
shape. When a source is located inside the cloak, it scatters waves as if it would be 
located some distance away from a small object. Moreover, the invisibility region 
actually hosts almost-trapped eigenstates. Mimetism is numerically shown to break 
down for the quantified energies associated with confined modes. If we further allow 
for non-isomorphic transformations, our approach leads to the design of quantum 
super-scatterers: a small size object surrounded by a quantum cloak described by a 
negative anisotropic heterogeneous effective mass and a negative spatially varying 
potential scatters matter waves like a larger nano-object of different shape. Potential 
applications might be for instance in quantum dots probing. The results within this 
paper as well as the corresponding derived constitutive tensors, are valid for cloaks 
with any arbitrary star shaped boundaries cross sections, although for numerical 
simulations, we use examples with piecewise linear or elliptic boundaries. 

1 Introduction 

Control of electromagnetic waves can be achieved through coordinate transformations 
which bring exotic material parameters [TJ [21 E] . Electromagnetic metamaterials within 
which negative refraction and focussing effects involving the near field can occur pEJ El [7J [8] 
can be understood in light of transformation optics [3]. 
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Recently, an electron focussing effect across a p-n junction in a Graphene film, that mimics 
the Pendry-Veselago lens in optics has been proposed [§]. The subsequent theoretical 
demonstration of 100 per cent transmission of cold rubidium atom through an array 
of sub de Broglie wavelength slits, brings the original continuous wave phenomenon in 
contact with the quantum world |10j . 

Other types of waves such as water waves can be controlled in a similar way using trans- 
formation acoustics [UJ, leading to invisibility cloaks for pressure waves thanks to the 
design of two-dimensional [T2l IT3] and three-dimensional cloaks [MJ [15] . It has been fur- 
ther demonstrated that broadband cloaking of surface water waves can be achieved with 
a structured cloak [16] . Interestingly, cloaking can be further extended to in-plane elastic 
waves [UJ [17] and bending waves in thin-plates [18] . 

In this paper, we focus our analysis on cloaking of quantum waves which involves spatially 
varying potentials and anisotropic effective mass of particles, as first proposed by the team 
of Zhang [19] and further mathematically studied by Greenleaf et al. [20]. We build up 
on the former proposal to render a quantum object smaller, larger, or even change its 
shape. Our point here is to apply the versatile tool of transformation physics in an 
area where the size of the object might have some dramatic changes in the physics: for 
instance, a quantum super-scatterer might enhance the interactions of quantum dots with 
the mesoscopic scale, thereby enabling quantum effects in metamaterials. 



2 Transformed governing equations for matter waves 

Following the proposal by Zhang et al. [19j, we consider electrons in a crystal with slowly 
varying composition: V = Et + U is the spatially varying potential, Ef, the energy of the 
local band edge and U a slowly varying external potential. In cylindrical coordinates with 
z invariance, and letting the mass density mo be isotropic diagonal in these coordinates, 
the time independent Schrodinger equation takes the form 

- • (m X W) + = E^ . (1) 

Here, h is the Plank constant and $ is the wave function. Importantly, this equation is 
supplied with Neumann boundary conditions on the boundary of the object to be cloaked. 
Let us consider a map from a co-ordinate system {u,v,w} to the co-ordinate system 
{x, y, z] given by the transformation characterized by x(u, v, w), y(u, v, w) and z(u, v, w). 
This change of co-ordinates is characterized by the transformation of the differentials 
through the Jacobian: 

, with 3 XU = |^4- (2) 
o(u, V, W) 

On a geometric point of view, the matrix T = J T J / det ( J) is a representation of the metric 
tensor. The only thing to do in the transformed coordinates is to replace the effective 
mass (homogeneous and isotropic) and potential by equivalent ones. The effective mass 
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becomes heterogeneous and anisotropic, while the potential gets a new expression. Their 
properties are given by [T5] 

Wi = m T^ , V = E + T~ l (V - E) , (3) 

where T^ 1 stands for the upper diagonal part of the inverse of T and T zz is the third 
diagonal entry of T. 

The transformed equation associated with the quantum mechanical scattering problem 
(JlT) reads 

— pi 

- 2 v • W= X V ^ + = m ' ( 4 ) 

where importantly the energy E remains unchanged and the wave function ty(x,y) = 
exp(iy/E(xki + yk 2 )) + ^d(x,y) with \/k\ + k\ = 1. We note that satisfies the usual 
Sommerfeld radiation condition (also known as outgoing wave condition in the context of 
electromagnetic and acoustic waves) which ensures the existence and uniqueness of the 
solution to ([5]). 

It is indeed the potential V and the mass density tensor m/_ (e.g. involving ultra- 
cold atoms trapped in an optical lattice as proposed in [T9]) which play the role of the 
quantum cloak at a given energy E. However, there is a simple correspondence be- 
tween the Schrodinger equation and the Helmholtz equation, the energy E of the former 
being related to the wave frequency u of the latter via tu = \[E (up to the normaliza- 
tion c = \/2/-</h, with c the wavespeed in the background medium, say vacuum). The 
present analysis thus covers cloaking of acoustic and electromagnetic waves governed by 
a Hemholtz equation. Correspondences bridging the current analysis with a model of 
transverse electric waves in cylindrical metamatrials are mf_ < — > d_ on the one hand and 
\/V' — E < — > lo on the other hand. 

3 Mathematical setup: Generalized cloaks for star 
domains 

This section is dedicated to a mathematical model generalizing the blowup of a point 
to a transformation sending a domain to another, thus making the latter inherit the 
same electronic, electromagnetic or acoustic behavior as the former, depending upon the 
physical context. Although in this paper we will restrict ourselves to cloaking regions in 
the plane or 3D Euclidean space, the transformation we propose can be readily extended 
to any star domain in W 1 , that is, domains with a vantage point from which all points 
are within line-of- sight. In particular, the transformation still preserves all lines passing 
through that chosen fixed point. 

Here is a description of the transformation in layman terms, but this could be formalised 
mutatis mutandis in very abstract mathematical settings by working directly with the 
divergence-form PDE of electrostatics [21] or as the Laplace-Beltrami equation of an 
associated Riemannian metric |22j. For simplicity, let us consider bounded star domains 
Dj in M n , n = 2, 3 with piecewise smooth arbitrary boundaries dDj, all sharing the same 
chosen vantage point, j = 0, 1, 2. We suppose D 2 contains Di which in turn, contains D Q . 
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Typically, D\ is the domain to be made to mimic D . The new transformation will be the 
identity outside D 2 , that is, in M. n \ D 2 , but will send the hollow region D 2 \ D to the 
hollow region D 2 \ D%, in such a way that the boundary dD 2 of D 2 will stay point-wise 
fixed, while that of Dq will be mapped to dD\. The hollow region D 2 \ D\ is meant to be 
the model for the cloak, endowed with Neumann conditions on its inner boundary dDi, 
and in which any type of defect could be concealed, but will still have the same electronic 
(resp. electromagnetic or acoustic) response as the region D with a potential wall (resp. 
infinitely conducting boundary or rigid obstacle) of boundary 8Dq. In practise, we may 
divide the domains Dj into subdomains, the part of whose boundaries lying inside dDj is 
a smooth arbitrary hypersurface. 

However, in such an ideal cloaking, there is a dichotomy between generic values of 
the energy E (resp. wave frequency w), for which the wave function must vanish within 
the cloaked region D\ \Dq, and the discrete set of Neumann eigenvalues of D% \ Do, for 
which there exist trapped states: waves which are zero outside of D\ \ D and equal to a 
Neumann eigenfunction within Di\D . Such trapped modes have been discussed in [22] 
when D vanishes. 




Figure 1: Construction of a generalized non-singular cloak for mimetism. The transformation 
with inverse |6p shrinks the region bounded by the two surfaces So and S 2 into the region bounded 
by S\ and S 2 . The curvilinear metric inside the carpet (here, an orange flower) is described by 
the transformation matrix T, see (f?|)-(fi^[j. This is designed to play the double role of mimesis 
and cloaking: any types of quantum objects located within the region D\ bounded by the surface 
S\ will be invisible to an outer observer while the region itself still scatters matter waves like 
an object Dq bounded by So (here, a yellow star). In the limit of vanishing yellow region, the 
transformation matrix T becomes singular on S± (ordinary invisibility cloak). 

The transformation is constructed as follows. Consider a point x of D\ \ Do with x = 
(x 1 , x 2 , ...) relative to a system of coordinates centered at the chosen vantage point 0. The 
line passing through x and meets the boundaries dD , dD\, dD 2 at the unique points 
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Xq, Xi, X2, respectively. We actually need the inverse x' n> x of the transformation, in 
the coordinates system, it reads x % = x\ + «j (x H — x\), where «j = x \~ x \ . In the 3-space 
with coordinates (x l ,x 2 ,x 3 ) = (x,y,z) we can write this transformation as 

x = x + a (x' — xi), with a = ^ 2 ~^° 

y = y + P{y'- vi) with = ggf 1 (5) 

z = z + 7 0' - zi), with 7 = a=a 

The cases of interest in this paper can all be considered as cylinders over some plane 
curves (triangular, square, elliptic, sun flower-like cylinders, etc.) Thus, we consider the 
transformation mapping the region enclosed between the cylinders So and S2 into the 
space between Si and S 2 as in Figure [TJ whose inverse is 

= x + a (x' — xi), with a = ^ 2 ~^° 

--yo + Hv'-yx) with /3 = g% xi (6) 

--z' 

The matrix representation of the tensor T _1 is thus given by 
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where the coefficients be expressed as 
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and finally the partial derivatives follows 

dx %2 — Xq %2 — x' dxo (X2 — xo)(x2 — x') dxi 

dx' x 2 — x\ x 2 — X\ dx' (x 2 — Xi) 2 dx' 

(x\ - x )(x' - xi) dx 2 
(x2 — xi) 2 dx' ' 
dx x 2 — x' dx (x 2 — x )(x2 — x') dxi 

dy' X2 — x\ dy' [x 2 — x\) 2 dy' 

(xi - x )(x' - xx) dx 2 
(x 2 - xx) 2 dy' ' 
dy y2-y'dy (y 2 - yo){y2 - y') dyi 



dx' y2 — yi dx' (y2 — yi) 2 dx' 

(yi - yo)(y' - yi) dy 2 



(z/2 - yi) 2 dx' ' 
dy _ z/2 - yo z/2 - y' dy (y2-yo)(y2-y')9y 1 



dy' 2/2 - 2/i 2/2 - 2/i dy' (y 2 - 2/i) 2 dy' 

(2/i - 2/o ) (2/ -2/1 ) dy 2 



(11) 



(12) 



(13) 



(14) 



(2/2 - 2/i) 2 dy' ' 

Now after having derived the general formulas for mimesis, we turn to the numerical 
simulations. From formulas ([6])- (14), in order to construct our cloak, we only need to know 
xo, xi, X2 and their respective derivatives. The explicit illustrations we have supplied to 
exemplify the work within this paper have boundaries whose horizontal plane sections 
are parts of ellipses (sunflower-like petal, cross, circle) or lines (parallelogram, hexagram, 
triangle). 



4 Mimetism for non-singular cloaks 

In Section [3] we presented the theoretical study of the mathematical model underlying 
our proposal for cloaks with any arbitrary star shaped boundaries cross sections, that 
perform mimetism as well as allowing invisibility. In the present section, we illustrate 
this by examples with piecewise linear or elliptic boundaries and provide their numerical 
validation. 



4.1 Formulas for piecewise linear boundaries 

If a piece of the boundary of a star domain Di is part of a line of the form y = a^x+hi, then 
clearly the line through the origin and a point (x', y') intersects this piece of boundary at 

to, yi) = ( , b%x , hlV , ) (15) 

2/ — QiX' y' — diX' 



and hence 

dxi biy' dxi bix' dyi a^hy' dyi ctibix' 



dx' (y' — ciix') 2 ' dy' (y' — a i x') 2, dx' (y' — aiX 1 ) 2 ^ dy' (y' — diX'Y 



(16) 
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Of course, in the case where this piece of boundary is a vertical segment with equation 
x = c, then the above intersection is at (c, c^). 

4.2 Formulae for piecewise elliptic boundaries 

We suppose here that a piece of at least one of our boundary curves Si is part of a 
nontrivial ellipse with equation of the form (x — a) 2 /c 2 + (y — b) 2 /d 2 where (a, b) is our 
ventage point. Of course a line (x(t),y(t)) = (a,b)+t(x' — a,y' — b) passing through (a, b) 
and a different point (x',y') intersects £j at two distinct points. For the construction, 
we need the point (xi,yi) of that intersection which is nearer to (x',y ; ) in the sense that 
Xi — a and y% — b have the same sign as x' — a and y' — b, respectively. 
We have 



Xi = a + 



x' 



VV - a) 2 /c 2 + (</' - WM 



b+ - y '~ b (17) 

v/(*' - a)7c? + (y' - WM 



This implies 



dxi 1 fx' — a) 2 



dx> ^(^ + (^ ^(^ + (^xf r2 ' 
(x' — a)(y' — 6) 

dy' 



c 2 + d 2 J "i 

dyi (x'-a)(y'-b) 



: 2 d 2 ' C « 



% I (y'-bf 



Or' / "'-"' 2 + fay ^ (^-a) 2 + fay 2 ^ i rf2 



(19) 
(20) 
(21) 



We can then apply formulae (|6|-(p!4|) to build any generalized cloaks involving boundaries 
of elliptic types, where the center of the ellipse is the ventage point (a,b). 

This has been used in Figure [2] (A), Figure [2] (C), Figures [5] (E) and Figure [9] (B). We 
note that outside the cloaks in Figure [2] (A) and [2] (C), the scattered field is exactly the 
same as that of the small disc of radius r = 0.195. When the radius r of the disc (that 



is, when c = d = r in equation (17)) tends to zero, the cloaks become singular and the 



plane matter wave goes unperturbed (invisibility). 
4.3 Squaring the circle 

In this section, we make a circle have the same (electromagnetic) signature as a virtual 
small square lying inside its enclosed region and sharing the same centre. In particular, 
its appearance to an observer will look like that of a square. As above, the transformation 
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(C) 




Figure 2: (Color online) A cross-like (A), a sunflower-like (C) and a hexagram (D) all mim- 
icking a circular cylinder (B) of small radius ro = 0.195, that is cq = do = 0.195 in equa- 
tion 111). The inner and outer boundaries of the petals are respectively parts of the ellipses 
x 2 /0.7 T + y 2 /0.2 2 = 1 and x 2 /0.9 2 + y 2 /0A 2 = 1 rotated by angle or § in (A) or 0, f , § or 
37r/4 in (C). The hexagram (D) is generated from an equilateral triangle with side 1.2. The en- 
ergy corresponding to the plane matter wave incident from the top is E = y/uj = yJlixcjX = 4.58, 
where A = 0.3 is the wavelength of a transverse electromagnetic wave in the optics setting with 
c the celerity of light in vacuum, normalized here to 1. To enhance the scattering, aflat mirror 
is located under each quantum cloak and obstacle. 



will map the region enclosed between the small square and the outer circle into the circular 
annulus bounded by the inner (cloaking surface) and outer circles, where the sides of the 
square are mapped to the inner circle and the outer circle stays fixed point wise. To do 
so, we again use the diagonals of the square to part those regions into sectors. Indeed, 
the diagonals provide a natural triangulation by dividing the region inclosed inside the 
square into four sectors. The natural continuation of such a triangulation gives the needed 
one. In each sector, we apply the same formulas as above, where (xo,yo) are obtained 
from the small square as in Section 4.1 and for both (x%, yi)and (x^, 2/2), we use the same 



For the numerical simulation, we use a small square 
0.2 and two circles radii R± — c\ — d\ — 0.2, R 2 = = d<i = 0.4 all centred 



formulas for ellipses in Section |4.2 
of side L 

at (a, b) = (0,0). So in both the uppermost and lowermost sectors (see Figure 3 (E) 
formulas (15) for the square become (xo,yo 

rightmost sectors, we have (xo,yo) : 

Riy'(x' 2 



V*^ij Vi 



{Rix'ix 12 + y 12 ) 



{Lo, 

y' 2 ) 



(Lo^7,Lo), whereas in the leftmost and 
For all sectors, formulas (17) now read 
i = 1,2. 
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The inner boundary of the cloak corresponds to a potential wall (with transformed Neu- 
mann boundary conditions, which also hold for infinite conducting or rigid obstacles 
depending upon the physical context). We report these results in Figure [2] and Figure [ij 
Some Neumann boundary conditions are set on the ground plane, the inner boundary of 
the carpet and the rigid obstacle. 

4.4 Star shaped cloaks 

In this section, we call star shaped a region bounded by a star polygon, such as a pen- 
tagram, a hexagram, a decagram, as in [41] . Star shaped regions are particular cases 
of star regions. The design of star shaped cloaks requires an adapted triangulation of the 
corresponding region. That is, a triangulation that takes into account the singularities at 
the vertices of the boundary of the region. So that, each vertex belongs to the edge of 
some triangle. To the resulting triangles, one applies the same maps as in Section [3] See 
e.g. Figure [2] (D). 

4.5 Finite Element Analysis of the cloak properties 

We now turn to specific numerical examples in order to illustrate the efficiency and fea- 
sibility of the cloaks we design. 

4.5.1 Comparison of backward and forward scattering of isomorphic cloaks 

Let us start with a comparison of both backward and forward scattering for the cloaks 
shown in Figure [2] (D). We report in Figure [3] the amplitude of the matter wave above 
and below the scatter (cloak and/or obstacle) along the x-axis respectively for y — 1 
and y = —1.2. We note the slight discrepancy between the curves, which is a genuine 
numerical artefact: we have checked that the finer the mesh of the computational domain, 
the smaller the discrepancy (which is a good test for the convergence of the numerical 
package COMSOL Multiphysics). The mesh needs actually be further refined within 
the heterogeneous anisotropic cloak and the perfectly matched layers compared with the 
remaining part of the computational domain which is filled with isotropic homogeneous 
material. We note that the yellow curve (corresponding to the pentagram, see Figure [2] 
D) is most shifted with respect to the black curve (corresponding to the small obstacle 
on its own i.e. the benchmark, see Figure [2] B). This can be attributed to the irregular 
boundary of the cloak as analysed in the case of singular star-shaped cloaks in [41J: we 
considered around 70000 elements for the mesh in all four computational cases reported in 
Figure [3] in order to exemplify the numerical inaccuracies. For computations with 100000 
elements, the yellow curve is shifted downwards and is nearly superimposed with the 
black curve. Moreover, the strong asymmetry of the yellow curve in the upper panel of 
Figure [3] vanishes for 100000 elements. We attribute this numerical effect to the artificial 
anisotropy induced by the triangular finite element mesh of the pentagram. 
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4.5.2 Analysis of the material parameters within a cloak 

Another interesting point in this paper, is that the metamaterial in our model is non- 
singular, therefore enabling the implementation of potentially broadband cloaks over a 
wide range of wavelengths. Moreover, the cloaked objects display a mimesis phenomenon 
in that, they are designed to acquire any desired quantum or electromagnetic signature. 

The result of the numerical exploration of the three eigenvalues A" 1 (a;, y) of the mate- 
rial tensor T _1 , of Figure [4] (E), is reported in Figure [5] First, using Maple software, we 
have purposedly represented those eigenvalues in a wider domain, including the cloaked 
region, despite the boundary conditions on the inner boundary of cloak. One clearly sees 
that all of the three eigenvalues take finite values away from even inside the cloak itself. 
The material tensor T -1 is hence nonsingular. Due to the fourfold symmetric geometry 
of the cloak (all sectors are obtained by a rotation of any fixed one), we only need to look 
at those eigenvalues in one sector. We note that 4 > A^ 1 > A^ 1 > 0.2 in accordance with 
the fact that the cloak should display a strong anisotropy in the azimuthal direction in 
order to detour the wave. For a singular cloak A^ 1 tends to zero on the inner boundary 
of the cloak, while A^ 1 tends to infinity. 

We also represent the finite element simulation of A^ 1 , which exemplifies the fourfold 
symmetry of the isovalues within the circular cloak, a fact reminiscent of the fourfold 
symmetry of the square which the cloak mimics. 

5 Generalized mirage effect and almost trapped states 

It is known that a point source located inside the coating of a singular cloak (i.e. a cloak 
such that xo = yo = in ^ leads to a mirage effect whereby it seems to radiate from a 
shifted location in accordance with this geometric transformation |23j . 

5.1 Shifted quantum dot inside the transformed space 

This prompts the question as to whether a similar effect can be observed in non-singular 
cloaks i.e. when Xq or y Q are nonconstant. As it turns out, the physics is now much richer: 
we can see in Figure [6] that when the source lies inside the coating, it only seems to radiate 
from a shifted location in accordance with ([6]), but it is moreover in presence of a small 
object of sidelength 2xq- This can be seen as a generalized mirage effect which opens 
many new possibilities in optical illusions. Indeed, Nicolet et al. have proposed to extend 
the concept of mirage effect for point sources located within the heterogeneous anisotropic 
coating of invisibility cloaks to finite size bodies which scatter waves like bodies shaped 
by the geometric transform [21]. This is in essence an alternative path to our proposal 
for mimetism. However, this mirage effect can be further generalized to non-singular for 
which a finite size body located inside the coating will now create the optical illusion of 
being another body in presence of some obstacle, which bears some resemblance with Fata 
Morgana, a mirage which comprises several inverted (upside down) and erect (right side 
up) images that are stacked on top of one another. Such a mirage occurs because rays of 
light are bent when they pass through air layers of different temperature. This creates the 
optical illusion of levitating castles over seas or lakes, as reported by a number of Italian 
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sailors, hence the name related to Morgana, a fairy central to the Arthurian legend able 
to make huge objects fly over her lake. 



5.2 Field confinement on resonances: Anamorphism fall down 

Another intriguing feature of singular cloaks is their potential for light confinement asso- 
ciated with almost trapped states which are eigenfields exponentially decreasing outside 
the invisibility region. Such modes were described in the context of quantum cloaks by 
Greenleaf, Kurylev, Lassas and Uhlmann in [20]. These researchers discovered that such 
modes are associated with energies for which the Dirichlet to Neunmann map is not de- 
fined i.e. on a discrete set of values. Here, we revisit their paper in light of non-singular 
quantum cloaks, that is when we consider the blowup of a small region instead of a point. 
Our findings reported in Figure [7] for a star shaped and a rabbit-like non-singular cloaks 
mimicking a small disc of radius 0.195 bridge the quantum mechanical spectral problems 
(panels (a) and (b)) to the scattering problems (panels (c) and (d)) in the following way: 
we first look for eigenvalues (i.e. quantified energies E) and associated eigenfunctions \I> 
of the equation Eq. Q in the class of square integrable functions on the whole space M 2 
(note that here, as the metric is non singular, there is no need to consider a weighted 
Sobolev space). Note however that the set continuity conditions on the inner boundary 
of the cloak, instead of Neumann ones. This provides us with a discrete set of complex 
eigenfrequencies, with a very small imaginary part (also known as leaky modes in the 
optical waveguide literature). We neglect this imaginary part and launch a plane wave on 
the non-singular cloak (whereby the invisibility region is also included within the compu- 
tational domain as we once again set continuity conditions on the inner boundary of the 
cloak) at the very frequency given by the spectral problem, see panels (c) and (d). We 
clearly see that the inside of the cloak hosts a quasi-localized eigenstate whose energy is 
mostly confined inside a star (panel C) and a rabbit (panel D), both of which actually 
scatter like a small disc. 

6 Generalized super-scattering for negatively refract- 
ing non-singular cloaks 

In this section, we take some freedom with the one-to-one feature of the previous trans- 
forms and allow for space folding. This means that \xq\ > \x2\ > \x\\ and \yo\ > 1 3/2. | > \yi\ 
in Q while the %i (resp. the yi) i=l,2,3 all have the same sign, thus making a and (3 
strictly negative real-valued functions. It has been known for a while that space folding 
allows for the design of perfect lenses, corners and checkerboards (121 1131 HU 12] ■ But it is 
only recently that researchers foresaw the very high-potential of space folding as applied 
to the design of super-scatterers (ISl 061 HU HHJ S9]. We generalize these concepts to 
mimetism via space folding. 

The mapping leading to the super-scatterer is shown in figure [8] We would like to empha- 
size that here we get not only a magnification of the scattering cross section of an object, 
in a way similar to what optical space folding does for a cylindrical perfect lens, but we 
can also importantly change the shape of the object. 
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We illustrate our proposal with a numerical simulation for a square obstacle surrounded 
by an anti-cloak in Figure |9](a), which mimics a larger square obstacle, see Figure ^b). We 
note that the large field amplitude on the anti-cloak's upper boundary can be attributed 
to a surface field arrising from the physical parameters with opposite signs on the cloak 
outer boundary (an anisotropic mass density in the context of quantum mechanics and 
an anisotropic permittivity in the context of optics). It is illuminating here to draw 
some correspondence with electromagnetic waves, as the anisotropic mass density (resp. 
permittivity) indeed takes opposite values when we cross the outer boundary of the cloak, 
and this ensures the existence of a surface matter wave (resp. a surface plasmon polariton) 
clearly responsible for the large field amplitude (in the transverse electric wave polarization 
i.e. for a magnetic field parallel orthogonal to the computational plane. We further show 
an example of a small circular obstacle mimicking a large square obstacle in Figure (9](a) 
and [9](b) . Once again, the large field amplitude on the outer cloak boundary comes from 
the complementary media inside and outside the cloak. We believe such types of mimetism 
might have tremendous applications in quantum dot probing, bringing the nano-world a 
step closer to metamaterials. 



7 Conclusion 

In this paper, we have proposed some models of generalized cloaks that create some 
illusion. We focussed here on the Schrodinger equation which is valid in a number of 
physical situations, such as matter waves in quantum optics. However, the results within 
this paper can easily be extended to the Helmholtz equation which governs the propagation 
of acoustic and electromagnetic waves at any frequency. One simply needs to insert the 
transformation matrix within the shear modulus and density of an elastic bulk (in the 
case of anti-plane shear waves), the density and compressional modulus of a fluid (in the 
case of pressure waves), or the permittivity and permeability of a medium (in the case of 
electromagnetic waves) [HI [12} [131 El 113 HE]. In this latter context, this means that in 
transverse electric polarization (whereby the magnetic field is parallel to the fiber axis), 
infinite conducting obstacles dressed with these cloaks display an electromagnetic response 
of other infinite conducting obstacles. In these cloaks, an electric wire could in fact be 
hiding a larger object near it. Actually, any object could mimic the signature of any 
other one. For instance, we design a cylindrical cloak so that a circular obstacle behaves 
like a square obstacle, thereby bringing about one of the oldest enigma of ancient time : 
squaring the circle! The ordinary singular cloaks then come as a particular case, whereby 
objects appear as an infinitely small infinite conducting object (of vanishing scattering 
cross-section) and hence become invisible. On the contrary, such generalized cloaks are 
described by non-singular permittivity and permeability, even at the cloak's inner surface. 
We have proposed and discussed some interesting applications in the context of quantum 
mechanics such as probing nano-objects. 

Obviously, one realizes that, when the inner virtual region D responsible for mimetism 
tends to zero, one recovers the case in [T] where the material properties are no longer 
bounded, as one of the eigenvalues of the mass density matrix tends to zero whilst another 
one recedes to infinity, as we approach the inner boundary of the coated region, see also 
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IZB. 

In this paper, we have played with optical illusions, trying to be as imaginative as 
possible in order to exhaust the possible geometric transforms we had at hand in Euclidean 
spaces (Non-Eulidean cloaking is a scope for more creative thinking [28]). It should be 
pointed out that while the emphasis of this paper was on quantum waves, corresponding 
non-singular cloaks in electromagnetism that have an inner boundary which is perfectly 
electric conducting (PEC) and scatter like a reshaped PEC object, were investigated in 
[HU EHl |56j EHl [57J |58j EH]. However these works focussed mostly on the reduction of 
the scattering cross section of a diffracting object, while the present paper explores the 
mimetism effect whereby an object scatters like another object of any other scattering 
cross section (and in particular reduced or enhanced ones). 

Metamaterials [5U] is a vast area with a variety of composites structured on the sub- 
wavelength scale in order to sculpt the electromagnetic wave trajectories, as experimen- 
tally demonstrated at microwave frequencies by a handful of research groups worldwide 
[5T| l52| 153] . Resonant elements within metamaterials are in essence man-made atoms 
allowing to mimic virtually any electromagnetic response we wish, and this is turn allows 
us to push the frontiers of photonics towards previously unforeseen areas. 
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Figure 3: (Color online) Upper panel: Profile of backward matter wave along the x-axis for 
y = 1 for a plane wave incident from the top as in Fig. ^for a cross-like (solid blue curve, see 
A), a sunflower-like (dashed red curve, see C), a hexagram (dotted yellow curve, see D) and a 
circular cylinder of small radius ro = 0.195 (solid black curve, see B); Lower panel: idem for 
profile of foreward matter wave for y = — 1; We note the large amplitude of the forward wave, 
due to the presence of a mirror below each nano scatterer at y = —1.2. The slight discrepancy 
between the curves is attributed to a numerical inaccuracy induced by the highly heterogeneous 
nature of the cloaks which is further enhanced by the irregularity of the boundary of the starshape 
cloak (dotted yellow curve, see D). The amplitude of the wave in both panels can be normalized 
to 1 (dividing throughout by 4.5 10~ 3 and by 5 10~ 3 in the upper and lower panels respectively) 
for comparison with Fig. [£| 
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Figure 4: (Color online) Left: A hollow parallelogram cylindrical region (A) scatters any 
incoming plane wave just like a much smaller solid cylinder (B) of the same nature. In (C) 
the same hollow parallelogram cylinder is designed to have the same response to waves as the 
small equilateral triangle (D) with side ^ip. Right, squaring the circle: metamaterials allow to 
make a circular cylindrical hollow region (inner and outer radii 0.2 and 0.4 respectively) and a 
small solid square cylinder (of side Lq = 0.2 and having the same center) equivalent, as regards 
their signatures and the way waves see them. In all these cases, the coated regions not only 
gain the same electromagnetic signature as any desired other object, but also serve as cloaks 
with nonsingular material properties, in fact the presence of any types of defects hidden inside 
them has no effect in the way they scatter waves. The energy corresponding to the plane matter 
wave incident from the top is E = ^JZj = y/2irc/\ = 4.58, where A = 0.3 is the wavelength of 
a transverse electromagnetic wave in the optics setting with c the celerity of light in vacuum, 
normalized here to 1. 
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Figure 5: ( A) -(C) are illustration of the graphs {(x,y, A^ )}, i = 1,2,3 of the three eigenvalues 
A^ 1 of the material tensor T _1 , in the uppermost sector of Fig. g] (E): (A) A^ 1 ; (B) A3 1 ; (C) 
A~2 . We note that, each of all those three surfaces are strictly above the plane z = even inside 
the cloak itself. Because all other sectors of the cloak are obtained by a rotation of the uppermost 
one, it suffices to study the eigenvalues o/T -1 in just one sector. Note that, the above surfaces 
were drawn using the MAPLE software. In (D) we represent the finite element computation of 
A 3 (see B) for all four sectors of Fig. g] (E) in the COMSOL MULTIPHYSICS package. We 
note the four-fold symmetry of the isovalues. 
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Figure 6: A non-singular cloak with two square boundaries of sidelengths 1 and 1.4 in presence 
of a quantum dot with energy E = 4.58 (resp. an electric current line source of wavelength 
A = 0.3 in the context of tranverse electric waves). (A-C) When the source is located a distance 
0.2 away from the cloak, it seems to emit as if it would be a distance 0.4 away from a square 
obstacle of sidelength 0.4. (b-d) When the source is located a distance 0.1 away from the inner 
boundary of the cloak (i.e. in the middle of the coating), it seems to emit as if it would be 
a distance 0.25 away from a square obstacle of sidelength 0.4, in accordance with |6p where 
xq = 0.2, x\ = 0.5 and X2 = 0.7. 
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(A) 




Figure 7: Left panel: Modulus of the fundamental eigenstates associated with quantified nor- 
malized energy E = y/oo = \f2nc/\ = 1.32 (resp. a transverse electric plane wave of wavelength 
A = 3.6 ^ for a non-singular cloak shaped as a star (A) and a rabbit (B) both of which mimic 
a small disc (of normalized radius 0.195, e.g. 195 nanometers) ; Right panel: Matter wave 
incident from the top on the quantum cloaks with a spatially varying potential V with compact 
support (i.e. vanishing outside the cloak) and energy E; The large amplitude of the field within 
the cloak is noted. 
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Figure 8: Construction of a generalized cloak with optical space folding for superscattering 
effect. The transformation with inverse |6p magnifies the region bounded by the two surfaces 
Sq and S% into the region bounded by Si and S2. We note that the transform is no longer an 
isomorphism. The curvilinear metric inside the carpet (here, an orange flower) is described by 
the transformation matrixT, see ffil-ftLfy. Any quantum object located within the region bounded 
by the surface S± scatters matter waves like a larger object bounded by Sq (here, a yellow star). 
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Figure 9: (A-C)Any obstacle surrounded by a square anti-cloak with square boundaries of 
sidelengths 1 and 1.4, scatters a plane matter wave of energy E = 4.58 (resp. a transverse 
electric plane wave of wavelength A = 0.3 in the context of optics) which is coming from the 
top like a larger square obstacle of sidelength 1.6. (b-d) Any obstacle surrounded by a circular 
anti-cloak with circular boundaries of radii 0.2 and 0.4 scatters a plane matter wave of energy 
E = 4.58 (resp. a transverse electric plane wave of wavelength A = 0.3J which is coming from the 
top like a larger square obstacle of sidelength 1. The large field amplitude on the upper boundary 
of the anti-cloak in (A) and (B) is noted. It can be attributed to some kind of surface matter 
wave (a surface plasmon polariton in the context of optics). 
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